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Nonlinear normal modes in electrodynamic systems: A nonperturbative approach 
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We consider electromagnetic nonlinear normal modes in cylindrical cavity resonators filled with a 
nonlinear nondispersive medium. The key feature of the analysis is that exact analytical solutions 
of the nonlinear field equations are employed to study the mode properties in detail. Based on such 
a nonperturbative approach, we rigorously prove that the total energy of free nonlinear oscillations 
in a distributive conservative system, such as that considered in our work, can exactly coincide 
with the sum of energies of the normal modes of the system. This fact implies that the energy 
orthogonality property, which has so far been known to take place only for linear oscillations and 
fields, can also be observed in the nonlinear oscillatory system. 
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I. INTRODUCTION 

The concept of normal modes is a milestone in the the¬ 
ory of linear oscillating systems and has had a significant 
impact on all fields of physics As is known, a lin¬ 

ear normal mode (LNM) is a free sinusoidal oscillation in 
a conservative dynamical system with constant parame¬ 
ters. In bounded distributed dynamical systems, an infi¬ 
nite but countable set of LNMs exists. In such systems, 
each LNM is characterized by its frequency and shape, 
and satisfies homogeneous partial differential equations 
(PDEs) of motion with given boundary conditions. A 
family of LNMs possesses the following important prop¬ 
erties, which allows one to solve a whole set of problems 
related to the calculation of free and forced motions in a 
linear system. 

1. Invariance. Each LNM can be excited indepen¬ 
dently of other LNMs by the specific choice of the initial 
conditions. 

2. Completeness. An arbitrary oscillatory process in 
the system can be expressed as a superposition of LNMs. 

3. Energy orthogonality. The total energy present in 
the system due to a free oscillatory process is the sum of 
the LNM energies. 

Since nature is nonlinear, LNMs can only be regarded 
as very useful mathematical models describing actual 
oscillations of nonlinear systems in the weak-amplitude 
limit. However, the following question naturally arises: 
Do strongly nonlinear systems admit such specific mo¬ 
tions that their properties allow one to consider them 
as nonlinear normal modes (NNMs), i.e., nonlinear gen¬ 
eralizations of the LNMs of the underlying linear sys¬ 
tems? An affirmative answer this question with respect 
to the lumped systems has been done in the seminal 
works of Lyapunov @ and Rosenberg [^-Q. Rosenberg 
defined an NNM as a vibration in unison of the mechan¬ 
ical system, i.e., a synchronous oscillation during which 
all the displacements of the material points of the system 
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reach their extreme values and pass through zero simul¬ 
taneously. A definition of the NNM in an autonomous 
distributed system in terms of the dynamics on a two- 
dimensional invariant manifold in phase space has been 
proposed by Shaw and Pierre Using the invariant 
manifold approach, they have also developed the tech¬ 
nique of asymptotic series expansions for constructing 
NNMs for a rather wide class of nonlinear 1 -|- ID au¬ 
tonomous systems. In the past decades, the concept of 
NNMs in mechanical systems has been studied exten¬ 
sively by a large number of workers (see, e.g., works [ 13 - 
[T3| and references therein). At the same time, NNMs in 
electrodynamic systems remain poorly studied. 

In this work, we present a nonperturbative approach 
to the concept of NNMs in an exactly integrable, non¬ 
linear 2 -b ID electrodynamic system. It should be em¬ 
phasized that all the forthcoming results are exact, i.e., 
no asymptotic expansions will be used. Our approach is 
based on the theory developed in Em- The nonlinear 
PDEs considered in [ia [isj l and herein depend explic¬ 
itly on the independent variable (radial coordinate) and 
can formally be regarded as 1 -I- ID nonautonomous sys¬ 
tems. Therefore, the approach proposed in Q, which is 
restricted to the autonomous systems, cannot be applied 
directly. In the present study, we define an NNM as fol¬ 
lows: The nonlinear normal mode (NNM) of a bounded 
distributed conservative system is an oscillatory motion 
in which all of the field quantities oscillate periodically 
in time with the same constant period in the whole vol¬ 
ume of the system. Each of the NNM fields must exactly 
satisfy the nonlinear PDEs of motion and the boundary 
conditions, and reduces to the corresponding LNM field 
in the weak-field limit. Note that a similar generalized 
definition of the NNM as an unnecessarily synchronous 
periodic motion in a mechanical system was proposed 

in El. 


In what follows, we construct the electromagnetic 
NNMs in cylindrical resonators filled with a nonlinear 
nondispersive medium and discuss their properties in de¬ 
tail. It will be shown that the considered NNMs, as 
their linear counterparts, exactly satisfy the first above- 
mentioned property of invariance and, what is quite re- 
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markable, the third (energy orthogonality) property. 

II. BASIC EQUATIONS 


The energy conservation law in the considered nonlin¬ 
ear medium can easily be derived from the field equa¬ 
tions o and ©. Multiplying o by and (|^ by H ^ 
after some algebra we obtain 


Consider electromagnetic fields in a bounded cylindri¬ 
cal cavity of radius a and height L. We assume that 
the z axis of a cylindrical coordinate system (r, 4>, z) 
is aligned with the cavity axis and limit ourselves to 
consideration of axisymmetric field oscillations, in which 
only the E^ and El^ components are nonzero. We will 
also assume that the cavity is filled with a nonlinear 
nondispersive medium in which the longitudinal compo¬ 
nent of the electric displacement can be represented as 
Dz — Dq + a“^eoei[exp(Q;i?2) — 1], where Dq, Ei, and 
a are certain constants. Since even powers of Ez are 
present in the series expansion of the medium does 
not possess a center of inversion. This is inherent in, 
e.g., uniaxial pyroelectric and ferroelectric crystals, pro¬ 
vided that the z axis is aligned with the crystallographic 
symmetry axis. With appropriately chosen constants Z?o, 
El, and a, such an exponential constitutive relation cor¬ 
rectly describes dielectric properties of actual media lack¬ 
ing a center of inversion in the case of weak nonlinearity 
where we can restrict ourselves to the quadratic (in Ez) 
correction term to the linear dependence of Dz on Ez 
(see [Hill for more details). Note that the model of 
a nonlinear distributed system presented herein and its 
generalizations were also discussed in works [l9l - l^ . In 
this case, the Maxwell equations read 


dw 


V-S = 0 


with the energy density 


(aE — 1) expiaE) + 1 77^ 

w = eoEi -5-h 


( 6 ) 

(7) 


and the Poynting vector 

S = -TqEH. (8) 

In the weak-field limit {\aE\ <C 1), expression (7) reduces 
to the well-known textbook formula w = t^£\E‘^j2 -f 

/ioi7V2. 

Based on formulas ®, the method proposed in |18l | 
makes it possible to easily construct exact periodic solu¬ 
tions of the system of equations m and starting from 
the LNMs which satisfy the linear wave equation ([5]). It 
is the purpose of the forthcoming analysis to demonstrate 
that the constructed solutions can actually be identified 
as NNMs of the considered electrodynamic systems. 


III. NONLINEAR NORMAL MODES IN A 
CYLINDRICAL CAVITY RESONATOR 


dH 

dr 


+ Ih = s{E)^, 


dE 


dH 


dr dt 


where E = Ez{r,t), H = H^{r,t), and 


£{E) = dDz/dE = eoEi exp(ai?). 


Analytical solution for the oscillations of the Eono type 

(1) in a circular cylindrical cavity with perfectly conducting 
walls and the nonlinear filling medium described by the 

(2) dynamic permittivity ([3]) has been found in [l^ and is 
given by 

A = A Jo(Knpe“^/^) cos(k„6»), 

(3) iZ =sin(K„d), (9) 


An exact solution to the system of equations ([T]) 
and m can be written in implicit form as [l8| 

£; = £ (pe“^/2^T + aZopiZ/(2e}/")) , 

H = Zo-ie}/" (p T + aZopiZ/(2£}/^)) . (4) 

Hereafter, p = r/a, r = f(eo£ipo)~^^^/«: and Zq = 
(po/co)^^^- The functions £ and EL describe the electro¬ 
magnetic field in a linear medium and satisfy the equa¬ 
tions 

d‘^£ 1 d£ d‘^£ 
dp"^ p dp dr"^ 

and 

d£ dn 

dp dr ' 


where A is an arbitrary amplitude factor, Jm is a Bessel 
function of the first kind of order m, Kn is the nth pos¬ 
itive root of the equation Jo{k) = 0, and 9 = t + 
aZopH/{2£^' ). The electric and magnetic fields are de¬ 
scribed by the implicit functions E{r, t) and H (r, t) which 
are solutions of system of two transcendental equa¬ 
tions. These implicit functions exactly satisfy Maxwell 
equations o and as well as the boundary conditions 

E{a,t)=0, |A(0,f)| < oo. (10) 

The initial conditions can be obtained by substituting 
T = 0 into formulas (O to give 

E = AJo(«„pe“^/2) 

at t = 0 and 


iZ(r,0) = 0. 


( 12 ) 
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For a sufficiently large index n such that n > n*{a), 
where n* is a certain integer, the functions E{r,t) and 
H{r,t) become ambiguous and solution ([9|), obtained 
without allowance for dispersion, becomes inapplica¬ 
ble [l^. For n < n* , implicit functions E and iJ 
given by formulas o describe continuous periodic os¬ 
cillations with the time period Tn = ‘I'KjiOn-, where 
Wn = Satisfy the NNM definition for¬ 

mulated above, and correspond to the NNMs of the i?onO 
(TMono) type in the cavity. Specifying the integer index 
n and imposing the initial conditions m and the 
motion of the system follows the exact solution © of 
the nonlinear boundary value problem for equations o 
and @, i.e., no oscillations with indices differing from 
n are excited. Hence, the considered NNMs satisfy the 
invariance property. 

Let us now consider some important features of NNMs, 
which were not pointed out in our previous work [l8l| . 
First of all, the electric field in these modes does not oscil¬ 
late in unison in the whole cavity volume, i.e., the ampli¬ 
tudes of the field at different spatial points can reach their 
extreme values and pass through zero at different instants 
of time. The same is true for the magnetic field. It is 
clearly seen in Fig. 3 in [l^ that there are no synchronous 
oscillations at different spatial points in the £^020 mode 
discussed therein. In this respect, the considered elec¬ 
tromagnetic NNMs differ from the well-known NNMs in 
lumped and I-flD distributed mechanical systems [1,0, 
and can be called “internally resonant” [11, It will be 
shown below that the degree of oscillation synchronism 
at different points of the cavity resonator depends on the 
shape of the cavity. 

We now calculate the total energy W stored in each 
NNM of the Eono type. To this end, one should substitute 
the implicit functions E and H given by formulas 
into © and integrate w over the cavity volume: 


L 27 r 1 


W = J j J w{p,T, a)pdpd(f)dz. (13) 


000 


A remarkable result is that the quantity W is indepen¬ 
dent of the nonlinearity parameter a and exactly coin¬ 
cides with the total energy of the corresponding 

linear Equo mode in the cavity resonator filled with a lin¬ 
ear medium that has the permittivity e = eo£i = const. 
The rigorous proof of this fact is given in the Appendix. 
Note that the quantity Wq is calculated analytically as 


i 

= ireQEia^L J + V?)pdp = 

0 

1 

= TTe^SiO^LA? COS^{hinT) J jQ{Knp)pdp + 

0 

1 

-\-sin^{Knr) J Jl{Knp)pdp = 


= ^coSia^LA^JliKn)- (14) 

It is worth also nothing that the fundamental frequencies 
ujn of the NNMs are independent of the field amplitude 
and the total energy, and coincide with the eigenfrequen- 
cies of the LNMs of the £onO type in the underlying linear 
system. In the next section, we will show that the above- 
mentioned notable features of the NNMs hold for an¬ 
other electromagnetic system described by equations o 
and (21). 


IV. NONLINEAR NORMAL MODES IN A 
COAXIAL RESONATOR 


Assume that a coaxial cylindrical inner conductor of 
radius 6 (0 < 6 < a) is inserted inside the cavity con¬ 
sidered in the previous section. The NNMs of the £onO 
type in the nonlinear coaxial resonator can readily be 
constructed using formulas (jd]) from the corresponding 
LNMs in the linear resonator with a constant permittiv¬ 
ity £ = eo£i (a = 0). The electric fields of the LNMs of 
the EonO type must satisfy equation ([S]) and the following 
boundary conditions on the perfectly conducting walls of 
the coaxial volume: 


£{b,t) = £{a,t) = 0. (15) 

The LNMs fields are given by 

£ — A [do P^^O (pn^ Jo (P'n )Fo iPnP^\ O.Os{^pnT ) , 

T-i = -A[Ji{pnp)Yt){pn) - Jo(/i„)yi(y:i„p)] sin(^„r),(16) 

where Ym is a Bessel function of the second kind of order 
m, and /i„ is the nth positive root of the equation 


Ml3p)Yo{p)-Mp)Yo{^p) = 0, (17) 


where (3 = h/a. 

Substituting functions (ITSl) into formulas o, we ob¬ 
tain an exact solution to the system of equations o 
and © in implicit form. It can be verified that the 
boundary conditions m remain valid for the implicit 
functions E{r, t) and i7(r, t) given by (0]) and (ITBl) . These 
implicit functions are periodic in time with the period 
Tn = 2TTlu}n, where a;„ = /i„(eo£i/io)~^^^a~^, since the 
transcendental equations (0]) are invariant with respect 
to the time shifts r —)• r -|- ITn with integer 1. Therefore, 
equations O and (1161) describe the fields of NNMs in the 
coaxial resonator with perfectly conducting walls and a 
nonlinear nondispersive filling medium. 

The fields of NNMs satisfy the initial conditions 


E = A 


MPnpe^^^^)Yoipn) - Jo{Pu)Yo{pnpe^^/^) 


(18) 

at t = 0 and H{r,0) = 0, and possess the invariance 
property. 

Let us turn to the results of some calculations by for¬ 
mulas o and (011). Figure 1(a) shows the snapshots of 
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FIG. 1. Electric and magnetic fields as functions of p (solid 
and dashed lines, respectively) in the n = 1 mode of the 
coaxial resonator with /I = 0.001 (a) and P = 0.999 (b) at 
times r = 0, ri = 27r/(7^i), T 2 — 27r/(5^i) and rs = 7r/(2/ii) 


the normalized electric and magnetic fields of the NNM 
of the i?oio type in the coaxial resonator with (3 = 0.001 
(n = 1 and /ii = 2.65...) as functions of p at fixed in¬ 
stants of time T. This case corresponds to a thin inner 
conductor (coaxial wire) inside the cylindrical cavity. For 
comparison, similar plots for the NNM of the Eqiq type in 
the cavity with /3 = 0.999 {pi = 3141.59...) are presented 
in Fig. 1(b). In the limit /3 —>■ 1, a coaxial geometry tends 
to a thin flat layer. Note that Fig. 1 corresponds to the 
case of strong nonlinearity where a A = 1. 

The presented plots show that the electric fields at dif¬ 
ferent spatial points do not oscillate in unison. However, 
it is seen in Fig. 1 that the degree of synchronism of the 
field in the NNM of the Flgio type inside the cavity turns 
out to be dependent on the value of /3. The oscillations 
are closer to synchronous ones for /3 = 0.001. 

The deviations of E and from their values corre¬ 
sponding to the Eqiq mode in a resonator with e = eo£i = 
const (a = 0) for /3 = 0.999 are more significant than 
those for /3 = 0.001, i.e., the nonlinear effects become 
more pronounced with increasing /3. In the limit /3 —1, 
the eigenvalue pn is close to an integer multiple of tt. 


This implies the more efficient interaction of harmonics 
of the eigenfrequency in the spectrum of each NNM. One 
may say that each NNM is more “internally resonant” in 
this case. 

The total energy of each NNM of the i?onO type is 
again independent of a and exactly coincides with the 
total energy of the corresponding LNM in the coaxial 
resonator filled with a linear medium (see the Appendix). 


V. MORE GENERAL OSCILLATIONS. 

ENERGY ORTHOGONALITY OF NONLINEAR 
NORMAL MODES 

In this section, we consider more general electromag¬ 
netic oscillations which correspond to the presence of an 
infinite set of NNMs in a cylindrical (noncoaxial) res¬ 
onator. Let us state the following initial and boundary 
conditions for the linear wave equation (5): 

d£ 

£:(p,0) = $(p), —(p,0) = 4'(p), 0<p< 1,(19) 

£’(1,t)=0, |f(0,r)|<oo, 0 < r < oo, (20) 

where $(p) and 41 (p) are given functions. The boundary 
value problem defined by ®, (HU), and (uni) describes free 
electromagnetic oscillations with the given initial field 
distribution in a cylindrical cavity specified by the rela¬ 
tions p = r!a < 1 and 0 < z < L, which is filled with a 
linear medium (a = 0). The solution to the linear bound¬ 
ary value problem specified by ©, (HU), and (uni can be 
found in a standard way by the method of separation of 
variables [l^. As a result, the functions £ and 'H are 
written as 


OO 

= X! ’^oinnp) [Bn cos(K„r) -I- C„ sin(K„T)] , 

n—1 

oo 

H(p,t) = - ^ Ji(k„p) [H„sin(«;„r)-C„cos(K„r)], (21) 

n—1 


where 


Bn = 


Jf{K„ 


p^{p)Jo{Knp)dp, 


Cn = 


^nJ\ (^n) 


1 

J p'i'{p)Jo{Knp)dp. ( 22 ) 


Substituting series m into formulas we obtain an 
exact solution to system of equations o and in im¬ 
plicit form. Such an implicit solution describes free elec¬ 
tromagnetic oscillations which correspond to the presence 
of an infinite set of NNMs in the nonlinear resonator. The 
implicit functions E{p, r) and iL(p, r) determined by for¬ 
mulas ® and (1^ satisfy the boundary conditions (1^ . 
but correspond to somewhat different initial conditions 
compared with m- 
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For example, let us specify the functions $ and 'h in 
the simple form 

^p) = A{l-p^), (23) 

vl-(p) = 0. (24) 

This leads to 

Bn = 8AK-%h{Kn)]-\ Cn = 0. (25) 

At the initial time r = 0, the electric field distribution 
E{p, 0) in the nonlinear resonator is defined by the tran¬ 
scendental equation 

OO 

A = ^5„Jo(K„pe“^/^), (26) 

n—1 

while the magnetic field iJ = 0 as in the “seeding” linear 
problem (see (1^ 1. 




T 

FIG. 2. (a) Initial distribution of the electric field E(p,0)/A 
(solid line) and the function ^{p)/A (dotted line), (b) Os¬ 
cillograms of the electric (E) and magnetic (H) fields of a 
nonlinear resonator at p = 0.75 (solid and dashed lines, re¬ 
spectively), calculated by formulas (4) and (21) for aA = 0.5 

The electric field distributions $(p)/A and E{p,0)/A 
in the linear and nonlinear cases, respectively, are pre¬ 
sented in Fig. 2(a). Figure 2(b) shows oscillograms of the 
field components E and E[ at the fixed point p = 0.75 for 
T > 0. Figure 2 was plotted for a A = 0.5. 

It is important to note that the exact solution obtained 
for such aA corresponds to single-valued continuous func¬ 
tions E and E[ (see Fig. 2(b)). While an NNM solution 


becomes ambiguous with increasing n for any fixed aA, 
the solution yielded by formulas (|1]) and m remains 
single-valued at moderate aA. This fact is stipulated by 
the rapidly decreasing series coefficients Bn- 

Despite the difference between the initial conditions 
given by (1^ and (1^ (see Fig. 2(a)), the energy for 
the field distribution (1^5)) in the linear cavity resonator 
{a = 0) coincides exactly with the energy of the field de¬ 
fined by (ESI in the nonlinear case where a 0 (see the 
Appendix). Due to the energy conservation, the same 
is true for an arbitrary time instant r > 0. This fact 
implies a quite remarkable energy orthogonality prop¬ 
erty of NNMs. The total energy W of the complicated 
oscillatory process, which is described by exact implicit 
solution o, with £ and % given by (1^ . is merely the 
sum of the NNM (or LNM) energies (fTTl) : 

CO 

W = A-^J2^nWt^ = 

n—1 

1 



0 


— -Q^oSiLa^A^. (27) 

The performed analysis shows that the energy orthog¬ 
onality property of the NNMs holds for a rather wide 
class of initial conditions (HID. One should only ensure 
the convergence of Fourier series (ED and the absence of 
ambiguity of the implicit solutions. 

For the electromagnetic fields governed by the non¬ 
linear system of equations (1) and (2), the principle of 
superposition does not hold and the NNM fields lack the 
usual orthogonality property which takes place for linear 
normal modes in resonators 0. Moreover, the interac¬ 
tion of the considered NNMs in the forced oscillations can 
result in complex nonlinear dynamics with a singular- 
continuous (fractal) Fourier spectrum [l^ . Therefore, 
the observed energy orthogonality property of the NNMs 
seems especially interesting. 

VI. CONCLUSIONS 

In this work, to the best of our knowledge, we have 
presented the first nonperturbative approach to the ba¬ 
sic properties of NNMs in a distributed nonlinear sys¬ 
tem. The approach does not require asymptotic expan¬ 
sions and provides a rigorous theoretical formulation of 
the NNM properties. It is to be emphasized that this 
formulation is not restricted to consideration of weakly 
nonlinear systems. 

In applying the developed approach, we have con¬ 
structed exact solutions to the electromagnetic fields of 
NNMs in cylindrical resonators filled with a nonlinear 
nondispersive medium. It has been shown that the field 
oscillations in the found NNMs are periodic in time, but 
are not synchronous at different spatial points. We have 
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established that the total energy of any NNM is indepen¬ 
dent of the nonlinearity parameter and exactly coincides 
with the energy of the corresponding LNM in the linear 
resonator. We have also obtained an exact solution which 
describes a more general oscillatory process correspond¬ 
ing to the presence of a countable set of NNMs in the 
nonlinear cylindrical resonator. Based on this solution, 
we have rigorously established the energy orthogonality 
property of the NNM helds. A very intriguing and phys¬ 
ically important issue, which naturally arises from the 
present analysis and still remains open, is whether en¬ 
ergy orthogonality of NNMs is a property inherent in the 
particular model of nonlinearity or can be extended to a 
wider class of nonlinear distributed systems. 


Consider for clarity a cylindrical (noncoaxial) res¬ 
onator. Making the change of variables and using (IA.4I) . 
from © and m one obtains the total energy 

W =2t: lO?La~^ x 
1 

X j{aE + 6-“-® - 1) X 
0 

x[l - aRSy2]RdR. (A.5) 

Here, we have also taken into account the boundary con¬ 
dition E = 0 at p = 1, implying that i? = 1 for p = 1. It 
is convenient to rewrite (IA.5I) as 
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W =2'Ke[)£ia^ L 




+—H ^(e —l)i?-|- 

la 


+ -ZR - 
a 2a ^ 


dR. 


(A.6) 


APPENDIX 

Let us show that the total energy of free nonlinear 
oscillations is independent of the nonlinearity parame¬ 
ter a and coincides with the total energy in the linear 
case (a = 0). Due to the energy conservation in a cav¬ 
ity with perfectly conducting walls and a nondispersive 
filling medium, it is sufficient to prove this fact for an 
arbitrary fixed time instant (say, r = 0). The electric 
field distribution E{p, 0) in the considered nonlinear os¬ 
cillations is defined by the transcendental equation 

E = Z{pe'^^/\Q), (A.l) 

while H{p,0) = 0. Introducing the notation R = 
pexp(aA/2), we have 

dR = {l + apE'p/2)e°‘^/^dp. (A.2) 

The derivative Ep can be found from (|A.1I) as 

E'p = [l- (A.3) 

Substituting (IA.3I) into (IA.2I) yields 

dp=[l- aRZ'j2]e-°‘^^‘^dR. (A.4) 


Integrating the first term in the integrand of (IA.6I) by 
parts and using the boundary condition f = 0 at i? = 1, 
we have 

1 1 

ZE'^R^dR = ^J Z^RdR. (A.7) 

0 0 

Integrating the second and third terms in the integrand 
of (IA.6I) by parts, one can find that the result of integra¬ 
tion cancels the last two terms in this integrand. Finally, 
we get 

1 

W = neoEio^L j Z^RdR, (A.8) 

0 

which coincides with the total energy of the linear oscil¬ 
lations. For a coaxial resonator, the proof is similar. 

In addition, it should be noted that the total energy 
of the radially localized field distributions vanishing for 
p —>■ 00 in an unbounded nonlinear medium, which can be 
the case for, e.g., cylindrical electromagnetic waves [l8l| . 
is also independent of the nonlinearity parameter a. 
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